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Abstract 

We computed the NLO corrections to twist-3, L — > T, flavor non-singlet amplitude 
in DVCS on a nucleon in the Wandzura-Wilczek approximation. Explicit calculation 
shows that factorization holds for NLO contribution to this amplitude, although the 
CN | structure of the factorized amplitude at the NLO is more complicated than in the 

■ leading-order formula. Next-to-leading order coefficient functions for matrix elements 

of twist-3 vector and axial-vector quark string operators and their LO evolution equa- 
tions are presented. 
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Introduction 

Deeply virtual Compton scattering (DVCS) [1, 2] on a nucleoli, 7*iV — > jN' , is perhaps 
the cleanest hard reaction sensitive to the skewed parton distributions (SPD). According to 
celebrated factorization theorems [3, 4, 5] the leading term in the 1/Q 2 expansion of the 
DVCS amplitude, where Q 2 is large virtuality of the hard photon, can be expressed in terms 
of twist- 2 skewed parton distributions. Twist-2 SPD's are related to matrix elements of 
non-local twist-2 quark and gluon string operators. These matrix elements contain complete 
information about nucleon structure seen by a highly- virtual electromagnetic probe. For 
that reason in recent years DVCS has been the subject of extensive theoretical investigations. 
First experimental data have also became recently available (see e.g. [6, 7, 8, 9]) and much 
more data are expected from JLAB, DESY, and CERN in the future. 

The amplitude of 7*V — > 7V scattering receives contribution from both DVCS and 
Bethe-Heitler processes. Extraction of DVCS from the experimental cross-section is there- 
fore not straightforward. The helicity analysis reveals that the interference term contains 
contributions of different photon helicity amplitudes weighted with sines and cosines of the 
angle between leptonic and hadronic scattering planes in the center-of-mass of the scattered 
photon and the nucleon [10, 11]. At present, the leading-twist DVCS amplitude has been 
studied to the NLO accuracy [12, 13, 14]. The analysis has revealed that the photon helicity is 
conserved by the tree-level contribution. Photon helicity-flip amplitude appears at the NLO 
level. The L — > T amplitude, corresponding to scattering of longitudinal virtual photon off 
nucleon, appears as twist-3, 1/Q correction. As the typical experimentally accessible values 
of Q 2 are by no means large, such leading power corrections may have significant effects on 
some of DVCS observables. In addition, twist-3 corrections typically scale as y/—i/Q, with 
t denoting the square of the momentum transfer, so the size of twist-3 corrections increases 
with t. As it follows, taking into account these corrections is mandatory for understanding 
continuation of the twist-2 part of the DVCS amplitude to t — 0. 

The LO contribution to the DVCS amplitude on the nucleon to the twist-3 accuracy has 
been calculated in [15, 16]. Twist-3 corrections to the DVCS amplitude arise in a twofold way. 
First, as usual for a hard exclusive process, the hard amplitude with an additional parton, the 
gluon, taking part in the hard collision is suppressed by one power of 1/Q. Such a genuine 
twist-3 contribution involves matrix element of three-parton operator in a nucleon state. 
Second, there is a so-called Wandzura-Wilczek (WW) contribution which arises from a hard 
configuration with minimal number of hard partons involved in hard collision. Formally, the 
WW contribution arises because operators with external derivatives w.r.t. total translation 
in a transverse direction give nonzero contribution in the DVCS kinematics. As the current 
phenomenology of power corrections is consistent with an assumption that matrix elements 
of three-parton operators in a nucleon are small [17], one can conjecture that the WW 
contribution can provide a rather accurate numerical description of twist-3 corrections. 

In this paper we have computed the NLO contribution to the WW twist-3, L — > T 
flavor non-singlet DVCS amplitude. Besides obvious phenomenological applications, there 
is a broader, theoretical interest in such a calculation. A factorization theorem for twist-3 
DVCS amplitude has neither been considered nor proven in the literature. As a consequence, 
although the direct calculation showed that the LO L — > T amplitude factorizes, there is no 
guarantee that factorization prevails in higher orders as well. At the same time there have 
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been no direct calculations supporting or disproving factorization of twits-3 contribution to 
DVCS beyond the leading order. The aim of this paper has been to provide an explicit 
example of NLO correction to twist-3 WW amplitude. Although computation of one-loop, 
flavor non-singlet amplitude may seem trivial, based on experience with NLO corrections to 
twist-2 amplitudes, calculation of NLO twist-3 contribution has turned out to be technically 
quite involved, warranting a more detailed discussion. Our calculation demonstrates that 
the amplitude does factorize at the NLO. It is tempting to interpret it as a hint that the 
factorization holds for this particular twist-3 DVCS helicity amplitude to all orders. 

Our paper is organized in the following way: in the first section we introduce our defi- 
nitions and notations. The next two sections are devoted to discussion of intricacies of the 
calculation: in the second section we derive convenient parametrization of matrix elements 
of vector and axial- vector quark string operators up to twist-3 accuracy in the WW approx- 
imation and in the third one we demonstrate technicalities of our approach by calculating 
LO twist-3 term in WW approximation. The fourth section contains discussion of our main 
result - the NLO correction to the WW L —> T amplitude. Finally, we conclude. 



DVCS amplitude on a nucleon 

Let p, p' and q, q' denote momenta of the initial and final nucleons and photons, respectively. 
The amplitude of the virtual Compton scattering process 

7*(?) + w(p)->7(g') + JV(p')> (1) 

is defined in terms of the nucleon matrix element of the T-product of two electromagnetic 
currents : 

T^ = ij d 4 x e -^'^(p'\T [J^{x/2)J^{-x/2)\ \p) , (2) 

where Lorentz indices fi and v correspond to the virtual, respectively the real photon. 

We shall consider the Bjorken limit, where — q 2 = Q 2 — > oo, 2{p • q) — > oo, with xb = 
Q 2 /2(p ■ q) constant, and t = (p — p 1 ) 2 <C Q 2 . We introduce two light-like vectors n, n* such 
that 



n • n — 0, n* • n* — 0, n • n* — 1. (3) 

We shall work in a reference frame where the average nucleon momenta P = \{p + p') and 
the virtual photon momentum q are collinear along z-axis and have opposite directions. Such 
a choice of the frame results in the following decomposition of the momenta [10]: 

— o 

P = n H n 

o 2 

q = -2£n + —n 

A = p - p = -2£n* + m 2 in + A± (4) 
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with m 2 = m 2 — 1/4, t = A 2 being the squared momentum transfer, and 



Q 2 - 1 

2i = 2 ^g 2 + 4ew • ( 5 ) 



Finally, £' is given by 



£' = = + 0(1/Q 2 ) , (6) 

+ t/Q 2 + + t/Q 2 ) 2 + 16f 2-a; B 



with ig = 7T— • 

We define the transverse metric and antisymmetric transverse epsilon tensors *: 

<% = <r-n*n* v -n v n*», e^e^nV^. (7) 

In the following, we shall use the shorthand notation for 

a + = a M n M , a" = a M n* M , (8) 

where a M is an arbitrary Lorentz vector. 

In this paper we focus our attention on the L^T DVCS amplitude where the virtual 
photon has longitudinal polarization. It appears first at the twist-3 level.* In order to 
calculate the leading contribution to this amplitude it is sufficient to expand (4) and (5) to 
twist-3 accuracy: 

P = ±(p + p>)=n*, A=p'-p=-2^P + A ± , 

q = -2£P + ^|n, q' = q-A = ^n-A ± . (9) 

The LO L^T amplitude has the form [15, 16] : 

( n _i_ AC p\v- f P u A a 1 1 r 1 r ~ ^ 

T o+ = ~ (Pq) + W ^J^dx {F a (x,OC + (x,0-iei p F^x,OC-(x,0} .(10) 



Note that the combination (g + 4£P)^ is proportional to the longitudinal polarization vector 
of the initial photon. In order to have exact gauge invariance of the amplitude we keep in 
(10) terms of order A ± /Q applying prescription g v " — > g v f + , p for the twist-3 terms in 

the amplitude. The skewed parton distribution F^x, £) and F^x, £) can be related to the 
twist-2 SPD's H,E,H : E with help of WW relations. We shall give explicit expressions for 
them in the next section. Term proportional to -jp^ is required to ensure electromagnetic 
gauge-invariance: the amplitude (10) is electromagnetically gauge invariant: 

qx: = u: q ' u = o (ii) 



'''The Levi-Civita tensor e^vap is denned as the totally antisymmetric tensor with 60123 = 1 
*We adopt here kinematical definition of twist i.e., terms suppressed by 1/Q are of twist-3. 
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formally to the accuracy 1/Q 2 - In the LO in QCD coupling coefficient functions are given 
by simple expression which have poles at x = ±£: 

c ± (x,t) = — l —— ± ' 



x — (, + ie x + !; — ie 

It was demonstrated in [18] that despite presence of these poles the convolution integral in 
(10) is well-defined and therefore factorization is not violated at the LO. 

One final comment is in order here. Note that using symmetry of the GPD's, see (33) 
below, the LO amplitude can be written in a shorter form as 

r °+ = i£ ^k + wl£ fa ;^(^ ( * ,{)+< ^ ( * ,{) )- (12) 



In the following we shall, for notational simplicity, always rewrite lengthy amplitudes in a 
similar way, making use of symmetries of corresponding GPD's. One can easily make use of 
these symmetries to transform expressions into the original form. 

Light-cone expansions of the matrix element 

Calculation of NLO corrections to the L^T amplitude requires knowledge of the nucleon 
matrix element of vector and axial-vector non-local quark string operators 

( P 'mx)r^(y)\p) (13) 

and 

(p'mxh^iyM (14) 

to the twist-3 accuracy but for arbitrary x and y. We shall restrict ourselves to the Wandzura- 
Wilczek approximation i.e., we neglect quark-gluon-quark operators arising in the twist 
expansion. As we shall show below, the result is given in terms of distribution functions 
which parametrize matrix elements of string operators restricted to the light-cone: 

r°° d\ 1 1 

F,(x,0 = j ^ -e-^p'm-XnhM-^lp) , (15) 

F,(x,0 = j ^ ^e-^(p'm-Xnh lll ^(--Xn)\p} . (16) 

To simplify notation, we have omitted in the above formulae the path-ordered exponential 
connecting the quark fields. 

As discussed in [16, 19], to the twist-3 accuracy the above matrix elements can be 
parametrized as: 

F^ w (x,0 = P^F(x,Z) + F^(x,£). (17) 
F™ w (x,0 = P„F(x,0+F ± „(x,0. (18) 
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Here we have explicitly separated longitudinal, twist-2: 



F(x,Z)=n'FY™(x,Z) = ((-Y + )KH + El)(x,Z)-((-))E(x,t), 



m 



(19) 



F(x,0 = n'F™^) = « 7 + 75»# (*,0 - «^R£(x,0 , 



(20) 



and transverse, twist- 3 components: 



Fl m (x,0 = -^^,0 + /_ 1 i ^(G>,o^ + (x,«,o + ^aG , >,0^-(^^0) • (22) 

A shorthand notation ((. . .)) denotes U(p') . . . U(p) and m is the nucleon mass. 

Functions and G^ 1 and the Wandzura-Wilczek kernels W±(x,u,£) have been intro- 
duced in Refs. [16, 19]. They are defined as : 



d 



d 



(23) 



(9m <9£. 



(24) 



±^(* > -O0(u >x)- 9(x < -O0(u < 
The flavor dependence in the amplitude can be easily restored by a substitution : 

q=u,d,s,... 

The light-cone expansion of the matrix elements (13,14) can be derived in a similar 
manner. To the twist-3 accuracy needed for the present calculation the result is 

(P + A/2\$(x)~f a ij(y)\P - A/2) = J* due i{Px){u -^- i(Py){u+0 {P a F(u, f) + F°(u, + 



(25) 



(26) 



iP* (x + y) p A p L F(u, + iP°tx - y) P G p 1 (u, 



(27) 



(P + A/2\^(x)^^{y)\P - A/2) = y 1 due i{Px){u -^- i{Py){u+ ^ {P a F(u,£) + F£(u,f)+ 



iP CT (z + y) P KF(u, + iP*(a; - y)„Gf(u, , (28) 



where now both x and y are off the light-cone. New distributions Gi(u,£) and Gi(u,£) can 
be expressed through the skewed distributions (17), (18) : 

G?(u,0 = uFZ(u,0-^efF ±k (u,0-l^l k ^±kF(u,0 (29) 

Gf(«,0 = «FL(«,0-e<%(«,0-^ A ±fc F(«,0 (30) 

Keeping in mind calculation of the L^T amplitude it is useful to recall symmetry properties 
of skewed parton distributions. Symmetry properties of twist-2 distributions H(u,£), E(u,£), 
H(u,£), E(u,£) with respect to interchange u <-> — u are given by: 

H(-u,£) = -H(u,£), E(-u,0 = -E(u,0, (31) 

H(-u, = H(u, 0, E(-u, = 0- (32) 

As it follows, 

F(-«,_0 = -F(«,0, F ±M (- M ,0 = -Fi>,0 (33) 

F(-u, = F(u, 0, = (34) 

In the following we shall often need Fourier transforms of the matrix elements (27) and (28). 
For example, for the vector matrix element one finds: 

/ Wf e ' i{k ' y) I (^I^ te )(P + A/2|^(y + ^)^(y-^)|P-A/2> = 

* du {p°F(u, + FZ(u, - P CT (f(u, + Gf («, 0^) } ^\2iP + fc') *W(tiP - k) 

(35) 

Similar result holds for the Fourier transform of the axial matrix element. 

Covariant calculation of the L — > T amplitude in the WW approximation 

The aim of this section is to present a method of calculation of the L — > T DVCS amplitude 
in the WW approximation which can be efficiently applied at the NLO. 

The crucial simplification arises due to the fact that according to (10) in the present 
kinematics the Lorentz indices \x and v of the L — > T amplitude point in the longitudinal, 
respectively transverse directions. As it follows, /i can be carried only by longitudinal vectors 
P and q and v has to be carried by transverse components of other vectors present in the 
problem. Symbolically, one can classify possible contributions as: 

(q or py x {A! , itfAl , Fl , G\{u, fl , FL GJ(«, } (36) 

Hence, it is necessary to compute only those terms which can be casted in the above form. 
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Figure 1: Graphical representation of the function H(x, y, z) 



The starting point for calculation is the standard definition of DVCS amplitude (2) which 
can be rewritten as: 



= if*xe^^ 2 <j/\T[J!^(x/2)JZ^(-x/2)]\p) = 
i J d 4 x e -*ta+9>/2 y d A yd 4 z(p'\^{y)H{x,y,z)^{z)\p) 



(37) 



Function H (x, y, z) is given by connected Feynman diagrams with amputated external fermion 
legs, see Fig. 1. Note that equation (37) is valid only in the WW approximation i.e. if twist-3 
contributions involving three points quark-gluon-quark matrix elements are neglected. Using 
Fiertz identity we can further transform (37) into: 

T^ = iJ d A x e-^'W 2 J d A y J d 4 z {H°(x,y,z){p'\iP(y)^(z)\p) 

-Ha(x,V, z){p'\4>(y)w5ip(z)\p) \ ■ 



Here we have introduced notation 



H V(A)( X ,V, Z ) = -7^\lAl<Tl5)H(x,y,z 



(38) 



(39) 



Using translational invariance and equations (27), (28) and (35) one easily obtains corre- 
sponding expression in momentum space. The final result for the twist-3 part of T^ v which 
gives rise to the L — > T amplitude reads: 

= J 1 du J d 4 k 8\uP -k) j d 4 k' <5 4 (2£P + k!) J d 4 k"d 4 (q + £P- k") 



{ \f 1M + W (<*(«,{) A + fma^ - 



■p 2 dk" 



H°(k,k', k") 
H%(k, k', k") \ . 



(40) 
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k"+k'/2 




k"-k'/2 



c-k'/2 



k+k72 



Figure 2: Momenta flow for Feynman diagrams of H(k,k', k 1 ' 



k"+k'/2 




k+k" 

k-k'/2 k+k72 




Figure 3: Leading order contribution to H(k, k', k") 

Functions H v ^(k, k', k") are Fourier transforms of functions H v ^(x, y, z): 

H v{A) (k, k', k") =i j d A x e~ lk " x J d 4 ye iy ( k+k 'M J dS e -«(*-*'/2) H v{A) (x, y, z) . (41) 

Resulting flow of momenta is shown in Fig. 2. Due to 5-functions in (40) one can put k = uP, 
k! = — 2£P and k" = q + £P after differentiation. Thus, calculation of the twist-3 DVCS 
amplitude in the WW approximation has been reduced to the calculation of hard parton 
diagrams and their derivatives with respect to external momenta. 

Above considerations are quite general and valid at any order of perturbation theory if 
the factorization is not broken. Let us now illustrate them with the example of calculation of 
the LO contribution to the twist-3 L — > T amplitude (10). In this case the hard amplitudes 
H v ^(k, k', k") are given by diagrams shown in Fig. 3. Simple calculation gives: 



H v{A) (k,k',k") 



^7^75)7^ + ^)7^ 
(k + k") 2 + i0 



(42) 



i.e. in the Born approximation H v , A Jk,k',k") does not depend on k'. Note that one need 
to evaluate only one diagram as contribution arising from the crossed one can be obtained 
using symmetry considerations and (33). Explicit calculation according to (40) yields: 

T o+ = £ du u+ *_. Q {±AlP»F(u, + [(« + Z)P" - <FK(u, + F"GT(u, 



+ -tel k A k ± P»F(u, - [(« + OP" + q"]ie1 k F k (u, + P^ ±k G k (u, 



(43) 
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+ crossed 



Figure 4: Feynmann diagrams for the NLO correction in WW approximation for H(k,k', k") 



Using the identity which follows from (29), (30): 

Gftu, + 1&1F(u, + <c^G*(«, - = (« - 

one can easily cast (43) into the form 
1 



(44) 



rp/lU 



2(Pg) 



(Pq) 



(45) 



Taking into account the contribution from the crossed diagram one finds that the term 
proportional to the integral f\duF"(u,£) cancels out and the final expression reads: 



F^u,0-ie v ±k Fl(u,0 



F^u^) + iel k Fl(u,0 



(46) 



which is precisely the expression given in (10). 

One comment is in order here. In the present calculation of the Wilson coefficient one 
does not need to apply QCD equations of motion (EOM) like it was done in Refs. [15] 
and [20]. EOM had been used at the previous stage to derive the light-cone expansion of 
the matrix elements (27) and (28). After it is done, the appropriate convolution of matrix 
elements with the hard amplitude according to (40) leads directly to the twist-3 L — > T 
Compton amplitude. 



NLO corrections to the WW twist-3 amplitude 

Technique described in the previous section can be applied without any modifications to 
the calculation of the NLO flavor non-singlet, quark contribution to the L — > T Compton 
amplitude. The corresponding Feynman diagrams shown in Fig. 4 are the same as for the 
NLO contributions to twist-2 amplitude. Nevertheless, besides trivial technical complications 
due to derivatives of corresponding diagrams w.r.t. external momenta, some new subtleties 
appear here which deserve a detailed discussion before final presentation of NLO results. 

A convenient method for regularizing loop integrals arising at the NLO is provided by 
dimensional regularization. We kept space-time dimension d = 4 — 2e and applied the 
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MS subtraction scheme. For the 75 matrix in D-dimension we used the t'Hooft-Veltman 
definition in terms of four antisymmetric gamma- matrices Powerful check of the resulting 
calculation is provided by the electromagnetic gauge-invariance of the final answer. We 
computed all possible contributions according to (36) and check that in the sum of all 
diagrams only such combination of momenta appears, g M + 4£P M , which is consistent with 
electromagnetic gauge-invariance. 

So far, it is a standard discussion familiar from calculations of higher-order coefficients 
in twist-2 amplitudes. New and interesting element appears when one considers a possible 
structure of the NLO answer. Obviously, proportionality of the L — > T amplitude in (10) to 
gM_l_4^pM i s fixed by gauge-invariance arguments and does not change when one goes from LO 
to NLO accuracy. On the other hand, the the LO amplitude depends only on combinations of 
F\_ and F*. Contributions involving F and F, which might in principle be present according 
to (40), cancel. This cancellation is, however, specific to the Born approximation and does 
not hold at higher orders. We therefore expect that the NLO expression for the L — > T 
amplitude has the form: 



rppuV 
J 0+ 



(g + 4£P)" 
2(Pq) 



97 + 



2tt 

C ± (x^;^ 2 /Q 2 ) = 



(Pq) j J-i 

AlC(x,t;iM 2 /<?)F(x,t) 
2 a s {fi 2 ^ 



dx \ C ± (x, £; fi 2 /Q 2 ) ( F ±a (x, + ie^Ft(x, 



iti p KC{x^^ 2 /Q 2 )F{ X) i) 



+ 



C F C { l\x^^ 2 /Q 2 ). 



2tt 



(47) 



Here we used symmetry properties of skewed parton distributions (33) in order to rewrite 
the contribution from crossed diagrams in terms of direct ones. 

Note now that term containing A°[F(x, £) may interfere with a corresponding piece of 
the L — > L amplitude Tj£ which appears at the NLO and describes the twist-2 transition 
between longitudinally polarized initial and final photons. In the massless QCD the corre- 
sponding contribution to the hadronic tensor Tqq has logarithmic singularity ~ In [q 2 /q' 2 ] as 
the virtuality of the final photon tends to zero. Of course, the L — > L scattering amplitude 
obtained by contraction of Tqq with polarization vectors of initial and final photons vanishes 
in this limit. In the Feynman gauge the contribution to Tqq arises from the box diagram 
[14]. From the general parametrization, valid for non-zero q' 2 : 



T ™ = JFqf ^ {Pq) ~ PV ) ( q ' U(Pq) ~ PV<1 ' 2 ) ALL (48) 

one finds that formally Tqq has a twist-3 part which can indeed contribute to T^ v with index 
fi transverse and v longitudinal. Retaining only the contribution singular as q' 2 goes to zero, 
one finds: 



rpfJLV 
1 f 



00 \tw-3 

A LL 



(q + A^PYi- A1)A LL , 



+ 



In 



In [q 2 /q' 2 } 

{i + x)2i ' (Z + X )(S- X ) 



2? 



(49) 



^Self-consistency of such requires to perform additional finite renormalization of the local axial current 
operator [21]. 
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Note that expression for All is finite as the space-time dimension d = 4 but instead it has 
singularity when q' 2 — > 0. We regulated it by a small non-zero value of q' 2 . In order to 
obtain the NLO L — > T amplitude the contribution (49) must be subtracted from result of 
calculation of diagrams in Fig. 4 in d = 4 — 2e dimensions and with q' 2 ^ 0. In practice, in 
the Feynman gauge only the box graph gives contribution singular as q' 2 — > 0. 

Explicit calculation according to the discussion above gives in d = 4 — 2e dimensions the 
following result for the coefficient functions C±\x, £; fi 2 /Q 2 ), C(x, £; fi 2 /Q 2 ) and £; ^ 2 /Q 2 ) 
in (47): 



C{x^-^ 2 /Q 2 ) 



x + i 

-5 - 2 In 



1 + 2 In 



2£ 



:i + e in[/i 2 /g 2 ]) 



2£ 



+ 4 In 



2£ 



j-x 
i + x 



In 



1 1 



;i + e in[/i 2 /g 2 ]) + 



2e 



2£ 



+ In 



Z + x 



2e 



+ 



2x 



In 



(50) 
2£ 



(51) 



C(x,Z;fi 2 /Q 2 ) 



1 1 



;i + 6in[/i 2 /g 2 ]) + 



e-x 2 



2 - In 



2£ 



2£ 



In 



2£ 



(52) 



Two important comments are of order here. First, in the above equations it is understood 
that £ has infinitesimally small imaginary part £ — > £ — «e. Moreover, as the combination 
Fj_ a + ie^ F^ 1S continuous at the point x = — £, the resulting integrals in (47) are well 
defined, in complete analogy with the situation in the LO amplitude. However, in order to 
complete proof of factorization of the L T amplitude at the NLO one has to demonstrate 
more, namely that the structure of singular terms is indeed such that they can be absorbed, 
following the standard procedure, into renormalization of the LO L — > T amplitude. 

In order to show this one has to derive one-loop evolution equations for the combination 
F^a + i^apFl which appears at the LO in (12). One possibility is to use definitions (21), (22) 
and the well-known evolution equations of the twist- 2 SPD's F and F see, for instance, [3]. 
Another way is to compute directly the one-loop evolution of the transverse components of 
the light-cone matrix elements (15) and (16). The result is: 



pL 2 -^(F ±a + iei p Ft)(x,{i) 



The evolution kernels are defined as: 



2vr 



-cv 



du 



V + (u,x,0(F ±a + iei p Ft)(u,0 



+U + (u,x,t)(AlF + ieLA' ± F)(u,t) 



V+{u,x,g) 



6>(£ < x < u) - 0(u < x < £) 



u — X 



+ 



6>(-£ < x < u) - 6{u < x < -£) 



u-x (£ + v,y 



(53) 

(54) 
(55) 
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0(£ < x < u) - 9(u < x < £) 
- 0(-£ < x < u) - 9(u < x < -£) 
Here the plus prescription is to be understood as 



(Z + x) 



{j + x){^ + u) 
4^ 2 (« + 2 



(56) 
(57) 



(X(u, x, £))+ = X(u, x,£)-6(x-u) J dzX(u, z, f) 



(58) 



As one can see the evolution equation of the WW contribution is not homogenous, i.e. the 
RHS of (53) receives contributions which are not given in terms of the combination of GPD's 
on the LHS. Precisely these terms generate contributions to the one-loop renormalization of 
the LO amplitude which exactly match singular terms arising from "bare" coefficients C and 
C. Indeed, with the help of equations (54) and (56) one finds that convolution integrals of 
the tree level coefficient function with kernels V+ and U + do reproduce the pole contributions 
in (50), (51) and (52) as it should be: 
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, (59) 
(60) 



Now, combining equations (50), (51), (52) with (59) one obtains the NLO Wilson coeffi- 
cients in their final form: 
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The above equations are our main result. They explicitly show that the NLO correction 
to the L — > T amplitude can be put into a factorizable form in a self-consistent way. 
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Conclusions 

In this paper we have demonstrated factorizability of the NLO correction to the twist- 
3, Wandzura-Wilczek, flavor non-singlet DVCS amplitude corresponding to scattering of 
longitudinal photon off nucleon. We have shown that, similarly to the LO, the singularity 
structure of the corresponding Wilson coefficients is such that the convolution integrals with 
SPD's are well defined. After the interference with the twist- 2 LL amplitude is properly 
taken into account all IR singular terms arising in the calculation of the L — > T amplitude 
can be absorbed into renormalization of the LO result in a usual manner. Besides the 
phenomenological applications of new Wilson coefficients found in this paper our result 
strongly suggests that factorization holds for the twist-3, L — > T DVCS amplitude beyond 
LO and possibly to all orders in the as expansion. 
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